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EFFECTS OF WING FLEXIBILITY AND VARIABLE AIR LIFT UPON WING BENDING MOMENTS 

DURING LANDING IMPACTS OF A SMALL SEAPLANE 1 

By Kenneth F. Merten and Edgar B. Beck 


SUMMARY 

A smooth-water-landing investigation teas conducted with a 
small seaplane to obtain experimental wing-bending-moment 
time histories together with time histories of the various pa- 
rameters necessary for the prediction of wing bending moments 
during hydrodynamic impact. The experimental results were 
compared with calculated results which include inertia-load 
effects and the effects of air-load variation during impact. The 
responses of the fundamental mode were calculated with the use 
of the measured hydrodynamic forcing functions. From these 
responses, the uring bending moments due to the hydrodynamic 
load were calculated according to the procedure given in R. & M. 
Xo. 2221. This comparison of the time histories of the experi- 
mental and calculated wing bending moments showed good 
agreement both in phase relationship of the oscillations and in 
numerical values. 

The effects of structural flexibility on the wing bending 
moment were large, the dynamic component of the total moment 
being as much as 97 percent of the static component. Changes 
in the wing bending moment due to the variation in air load dur- 
ing impact were of about the same magnitude as the static 
water-load component. 

INTRODUCTION 

Recent trends in the design of aircraft have led to an 
important increase of the stresses produced in wings by land- 
ing impacts. Two significant factors contributing to these 
increased stresses during landing are an increased proportion 
of the airplane weight in the wings and an increased struc- 
tural flexibility, since, in most cases, these factors have caused 
the ratios of the times to peak of the applied landing loads 
to the period of the fundamental mode to approach a critical 
value. 

Several simplified methods have been developed for deter- 
mining the inertia loads in wing structures during landing 
impacts, and studies have been made of the landing-impact 
inertia loads in simplified structures with the use of the 
principles of these methods. Although experimental investi- 
gations have been made to determine the magnitudes of 
inertia loads in actual airplane structures, little correlation 
of theory and experiment has been made concerning the 
nature and magnitude of inertia loads in airplane wings 
during actual landing impacts. 


Another aspect of the problem of wing loads during land- 
ing is the variation of air load due to changes of attitude 
and flight path during impact. The importance of this 
change in air load has been the subject of some speculation, 
but little investigation. 

In order to evaluate the importance of the various com- 
ponents of the load, including dynamic effects and variation 
in air load, data were obtained during full-scale landing tests 
of a small seaplane to provide a comparison of actual wing 
loads with those predicted by a simplified method 
(reference 1). 

The present report gives a comparison of the theoretical 
and experimental wing loads, in the form of time histories of 
the wing bending moments, and discusses the contributions 
of each of the components of the moment (static water-load 
moment, dynamic water-load moment, and air-load moment) 
to the total. The static and dynamic components of the 
total moment were calculated and combined according to 
the procedure of reference 1, with the responses of the funda- 
mental mode being calculated from the recorded time histories 
of the applied forces. The air-load component was calcu- 
lated by a simplified method which is described in this report. 

SYMBOLS 

Cl lift coefficient 

g acceleration due to gravity (32.2 ft/sec 1 ) 

L lift, pounds 

M bending moment in wing, pound-inches 
n load factor, multiples of g 
S wing surface area, square feet 

T duration of impact, seconds 

ti time to peak of applied load, seconds 
f, natural period of fundamental mode, seconds 
u dynamic response factor, ratio of maximum total 
water-load wing bending moment to maximum 
static water-load bending moment 
V velocity of seaplane, feet per second 

a angle of attack, degrees 

y flight-path angle, degrees 

A prefix denoting change 

p density, slugs per cubic foot 

I t trim angle, degrees 


t Supersedes N"ACA TN 3063, “A Comparison of Theoretical and Eiperimental Wing Bending Moments during Seaplane Landings’' by Sennetb F. Merten, Ja*< L. Rodrigues, 
and Edgar B. Bed;, 1950 


221 



222 


REPORT 1013 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Subscripts: 

c corrected for air load 
h horizontal 

n normal to keel 

0 at time 'of water contact 

p parallel to keel 

r recorded 

T total 

v vertical 

max maximum 

APPARATUS AND INSTRUMENTATION 

The airplane used in the present investigation was a small 
two-engine seaplane (fig. 1). Pertinent information about 
the seaplane is given in table I, and additional information 
may be obtained from reference 2. The frequency and shape 
of the fundamental wing bending mode were found from 
ground vibration tests and are given in table II and figure 2. 
The spanwise weight distribution is also given in table II. 

The trim variation was measured with a gyroscopic trim 
recorder mounted in the cabin floor. The airspeed was 
measured with an NACA airspeed recorder, pitot-static-tube 
type, mounted above the cabin. Accelerations of the center 
of gravity were obtained from an NACA optical-recording 
three-component accelerometer mounted securely in the 
fuselage near the center of gravity. The time of contact 
was determined from a water-contact indicator located on 
the keel at the main step. The hull immersions were deter- 
mined from pressure gages installed along the bottom of the 
hull. The wing bending moments were measured by means 
of a strain gage mounted on the wing main spar 9 inches 
from the center line of the seaplane (hereinafter referred to 
as station 9) . 

The estimated accuracies of the experimental data based 
on calibration, instrument, and reading error are as follows: 


Horizontal velocity, I'*, feet per second ±4 

Trim angle, r, degrees ±0. 25 

Load factor, n, multiples of g ± 0. 2 

Initial wing lift, Lo, multiples of g ±0. 05 

Total wing bending moment, M t, pound-inches ±0. 05X 10 s 


TEST PROCEDURE 

The landing-impact tests were made in smooth water. 
During these landing tests, airspeed, trim variation, center- 
of-gravity accelerations, and wing-spar bending moments 
were recorded. The landings were made at horizontal veloc- 
ities ranging from 95.4 to 112.0 feet per second, trim angles 
ranging from 3.00° to 7.83°, and initial flight-path angles 
ranging from 2.0° to 4.4°. The resulting maximum center- 
of-gravity accelerations normal to the keel line ranged from 
l.lOjf to 1.96 j 7 and the duration of the impacts varied from 
0.63 to 0.87 second. The times to peak of these normal 
accelerations ranged from 0.10 to 0.36 second. Values of 
these parameters and other pertinent information for all the 
tests are presented in table III. 



1 NACA op tlcal -record Eng three-oomponent 3 Pressure sages and water contacta 

accelerometer 

2 Gyroscopic trim recorder 4 Water-contact Indicator 

Fioi-bi l. — i Seaplane and instrumentation. 

Table I 

GENERAL INFORMATION ABOUT SEAPLANE USED IN 
LANDING-IMPACT INVESTIGATION 


Approximate flying weight during tests, lb 19, 200 

Stalling speed (flaps down), fps 94 

Wing area, sq ft 780. 6 

First natural frequency, ops 4. 76 

Second natural frequency, cps 13. 0 



FinusE 2 . — Experimental fundamental bendlng-mode shape of wing semlspan. Frequency, 
4.76 cycles per second. 


Table II 

WEIGHT DISTRIBUTION AND UNITIZED MODE SHAPE 
FOR FUNDAMENTAL BENDING MODE (4.76 cps) 

OF WING SEMISPAN 


Distance Irom 
center line 
(in.) 

Associated 

Btation 

weight 

(lb) 

Unitized 
mode shape 

a 


-0.045 

31 

881 

-.044 

75 

2,007 
fig 076 

-.026 

87.7 

—.023 

110 

SSI 

-.004 

170 

116 

.063 

210 

102' 

.110 

230 

88 

.190 

290 

181 

.270 

330 

64 

.370 

370 

53 

.490 

410 

43 

.626 

440 

18 

.730 

477.7 

40 

.860 

510 


1.000 

Total 8,600 
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Table III 

FLIGHT AND IMPACT VARIABLES 


Ron 

ffps) 

(fps) 

ttn) 

flW) 


Bn 

•mma 




T 

(sec) 

h 

tj* 

(lb- In.) 

u 

(deg) 

(deg) 

(deg) 

(deg) 

a t.„ 

(deg) 

I 

mu 

4.91 

5.1 

8.37 

1*06 

2.03 

0.264 

LOO 

0.506 

0.69 

3.72 

0.573X10* 

L47 

7.83 

L30 

10.66 

5.30 

4.7 

a 

95.4 

3.55 

3.5 

7. 29 

L.21 

L26 

.181 

L 19 

.259 

.87 

6.28 

.250 

L23 

5.01 

5.07 

8.11 

2.10 

4.2 

3 

100.5 

3.46 

7.2 

11.23 

L61 

L 69 

.209 

L50 

.340 

.65 

3.43 

.410 

L 53 

5.72 

5.21 

9.82 

4.36 

6.5 

4 

112.0 

3.83 

8.6 

13.54 

LfiO 

1. 94 

.255 

1.89 

.355 

.63 

L 90 

.483 

L97 

3.00 

9.31 

7.39 

4.45 

7.0 

5 

mo 

3.54 

4.4 

8.10 

LIS 

L24 

.183 

L 14 

.285 

.77 

4.36 

.345 

L 42 

6.07 

4.2 8 

8.45 

3.29 

4.4 

6 

104.1 

3.42 

3.6 

6.90 

L. 10 

L 19 

.208 

LQ7 

.310 

.80 

6.S5 

.335 

L 13 

7.39 

3.30 

9.34 

2.67 

3.8 


THEORY 

WILLIAMS’ METHOD 

In reference 1 a method was proposed for calculating the 
dynamic effect of an impulsive load applied to an elastic 
structure. Basically, the method follows classical normal- 
mode vibration theory by considering the total response 
of the elastic structure to a forcing function at any instant 
to be the summation of the responses of all of its normal 
modes at that instant. However, a unique feature of the 
method is that the total response of each mode is separated 
into a static and a dynamic component, and the stress due 
to the sum of the static components of the responses of all 
the modes is found in one calculation by rigid-body analysis. 
This stress is referred to as the static-load stress. The 
stress of each mode due to its dynamic component of re- 
sponse is found separately. The total stress is the sum of the 
static-load stress and the dynamic components of stress 
for the significant modes. Time histories of stress are found 
in these calculations and thus phase relationships of the 
modes are considered when the stresses for each mode are 
added. 

AIR-LOAD VARIATION 

In determining the effect of air-load variations on the 
wing bending moment during impact, the air load was 
assumed to change instantaneously with change in angle of 
attack and the rate of change in air load was assumed to be 
slow enough to neglect excitation of structural dynamic 
response. Also, the ratio of air-load bending moment to 
lift was assumed to be constant for all impacts. On the 
basis of these assumptions, the air-load moment at any time 
during the impact M was related to the lift L and the ratio 
of bending moment to lift at time of contact MJLo by the 
following equation: 



An average value of Mj/Lq, obtained from the measured 
values of lift and moment at time of contact, was used in 
this equation. The variation of lift L was determined 
from time histories of the angle of attack and velocity in 
conjunction with a lift curve (C L plotted against a) determined 
for the wing with flaps deflected. This lift curve was 
determined by fairing a line with the slope of the available 
lift curve without flaps deflected through the values of 
Cl computed for all the impacts from the measured velocity, 


angle of attack, and cenfcer-of-gravity acceleration at tune of 
contact. With the variation in. total air-load moment thus 
determined, the change in air-load wing bending moment 
during impact Ail / was found from the following equation: 

A 

CALCULATIONS AND RESULTS 

The variations of the wing angle of attack aud velocity 
during impact necessary for computing the changes in air-load 
bending moment were determined for each impact from the 
recorded data in the following manner. The accelerations 
normal and parallel to the keel line, obtained from the 
three-component center-of-gravity accelerometer, were plot- 
ted (fig. 3). After the trim-angle variation (fig. 4 (a)) was 
taken into account, these accelerations were resolved into 
vertical and horizontal components. Integration of the 
time histories of these accelerations over the duration of 
the impact produced time histories of the changes in vertical 
and horizontal velocities. Since the vertical velocity at the 
time of contact was not accurately known, the initial 



(a) Acceleration parallel to seaplane keel line during Impact. 

(b) Acceleration normal to seaplane keel line during Impact. 
Flora* 3.— Typical accelerometer records and bakings. 
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velocity was determined so that integration of the time- 
history curve over the duration of the impact resulted in a 
final vertical displacement of zero. (The duration of the 
impact is defined as the interval between the time of contact 
and the time when the center of gravity again reached its 
initial height above the mean water line. The instant of 
contact was found from a water-contact indicator on the 
step and the time history of the center-of-gravity displace- 
ment was determined from the times of immersion and 
emersion of the hull pressure gages, the fixed location of the 
center of gravity relative to the step, and trim-angle time 
history.) Integration of the time-history curve of the 
corrected vertical velocity from time of contact to the time 
of zero vertical velocity dete rmine d the maximum displace- 
ment of the center of gravity. The maximum displacements 
determined in this manner for all the impacts agreed within 
experimental error with the maximum displacements cal- 
culated from the hull pressure gages. With the use of the 
corrected vertical-velocity and horizontal-velocity time 
histories, time histories of the flight-path angle y and the 
resultant velocity were computed. From the time histories 
of trim angle r (fig. 4 (a)) and flight-path angle y (fig. 4 (b)), 
the time history of the angle of attack a was computed 
(fig. 4 (c)). 

With the use of the time histories of angle of attack and 
resultant velocity, the changes in bending moment in the 
wing at station 9 due to the changes in air load were deter- 
mined for each impact in the manner described in the 
section entitled “Air-Load Variation” and are presented in 
parts (a) of figures 5 to 10. 



(a) Trim-armie variation daring Impact 

(b) Flight-path angle variation during Impact. 

( 0 ) Angle-of-attack variation during Impact. 

Fwubk 4.— Typical time histories of pitching motion based on recorded and oalculated data. 



(a) Components of wlng-bending-moment changes during Impact. 

(b) Total wing bending moment during Impact. 

Fiot'B* i.— Wlng-bendlngHnoment time histories daring Impact; ran I. 



i 

£ 


I' 

-5 


a / cy- 

— 

— — 

— 

,--L 

— T7 

evel-fl 

ighf mi 

vnent 

— 












,.T 

of a! th 

eorefk 

:a! mor 

lent 


s' 


A 



S' 








/ 



Vc 

fa! exf. 

lenmer 

ita! me 

menf 












(b) 










./ 


.5 


.3 


.7 


.8 


.4 .5 .6 

Time, sec 

fa) Components of wlng-bending-moment ohangej daring Impact. 

(b) Total wing bendiDg moment during impact. 

FiauEB 5.— Wlng-bending-moment time histories during Impact; run 2. 
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(a) Components of wing-bendlng-moment obanges daring Impact. (a) Components of wing-bendlng-moment changes daring Impact. 

(b) Total wing bending moment during impact. (b) Total wing bending moment dnrfng impact. 

Flora* 7. — Wing-bendlng-moment time histories during Impact; run B. Flora* 9. — WTng-h pndinp - rrmmpn t time histories 'hirin g- impact; run A 



(a) Components of wing-bendlng-moment cbanges daring Impact. 

(b) Total wing bending moment during Impact. 

Flora* A— Wing-bendlng-moment time histories during Impact; run 4. 
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(a) Components of wing-bendlng-moment changes during Impact. 

(b) Total wing bending moment during impact. 

Flora* 10.— Wlng-bending-moment time histories during Impact; run 8. 
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The procedure of reference 1 was used to compute the 
bending moments because it provides a convenient means of 
applying the principles essential to a dynamic-loads analysis 
which results in time histories of the wing bending moments. 
The forcing function for each impact was determined from 
the nor mal acceleration measured in the hull by an accelerom- 
eter located near the seaplane center of gravity. Because 
of the manner of connection between the wing and hull 
(fig. 1), the measured acceleration was not appreciably 
affected by the oscillations of the whig. By including the 
effect of the varying wing lift on the center-of-gravity 
acceleration, the acceleration normal to the keel line due 
to tile hydrodynamic force only w T as determined. From this 
acceleration and the mass of the seaplane, the hydrodynamic 
forcing function was calculated. The dynamic responses 
of the significant modes to this forcing function were computed 
by a recurrence method developed at the Langley Structures 
Research Division in which the actual forcing function 
is used without approximation. Only the dynamic effects 
of the fundamental bending mode were included in the 
final results because calculations showed the dynamic effects 
of the second symmetrical bending mode to be negligible. 
This observation was borne out by the absence of higher 
mode effects on the strain-gage records. The calculated 
time histories of the static water-load and dynamic water- 
load components of wing bending moment at station 9 are 
presented in parts (a) of figures 5 to 10. The spanwise 
bending-moment distribution for the fundamental-mode 1 g 
inertia loading calculated as set forth in reference 1, a 
1 g static water loading, and a level-flight loading are plotted 
in figure 11. The values of the bending moment at station 
9 used in the application of the method of reference 1 wore 
obtained from this plot. 

These three components of wing-bending-moment changes 
obtained in this manner for station 9 were combined and 
added to the wing bending moment existing at the instant 
of contact. This total theoretical wing bending moment is 
presented in parts (b) of figures 5 to 10 together with the 
wing-bending-moment variation measured by the strain 
gage at station 9. 

DISCUSSION 

Comparisons between the total theoretical and experi- 
mental wing bending moments are presented in parts (b) 
of figures 5 to 10. The comparisons are made only for wing 
station 9 because the bending moments in the outer wing 
section were so small as to be of the same order as the 
estimated error. Only the dynamic effects of the funda- 
mental bending mode were included, in the final results 
because calculations showed the dynamic effects of the 
second symmetrical bending mode to be negligible. This 
negligibility was borne out by the absence of higher-mode 
effects on the strain-gage records. The comparisons show' 
the predicted values to be in good agreement with the experi- 
mental values. As can be seen from the figures, the phase 
relationships between the theoretical and experimental 


1.4x10 4 



Fjgvbe II.— Level-flight and landing-impact bending moments. 


values are consistently good, and the maximum changes 
from initial conditions show a range of error of 5 to 28 per- 
cent based on the experimental values of the maximum 
changes in wing moment. These results indicate that when 
the three components of moment are included in the theory 
good agreement is obtained. 

A comparison of the level-flight bending moment and the 
change in bending moment due to the heaviest impact of 
the tests shows that the maximum change in bending mo- 
ment accompanying a downward motion of the w’ing w’as 
approximately 50 percent of the level-flight moment (figs. 1 1 
and 5 (b)). This maximum bending-moment change was 
produced by a 2.030 impact. When any differences in 
the response factor and any change in air-load bending 
moment are neglected, landings of over 4 g would be required 
for the downward motion of the w'ing merely to start stress 
reversal in the wing. Similarly, even if the maximum 
dynamic w r ater-Ioad bending moment (fig. 8 (a)) caused by 
a 1.940 impact were to be twice as large for a 4 g impact, 
were to exist after the water load was removed, and were to 
be superposed on a level-flight moment, the moment pro- 
duced by the upward motion of the wing w r ould still be less 
than twice the level-flight moment. Therefore, the change in 
bending moment due to a landing impact is unimportant 
in this seaplane, insofar as this change will not produce 
critical stresses at the wing root. This unimportance may 
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be largely attributed to tbe fact that the fundamental 
mode 1 g inertia loading is relatively small as compared 
with a level-flight loading (fig. 11). 

The effects of structural flexibility on the computed bend- 
ing moments can be seen by comparing the static and dy- 
namic components of the water-load bending moment in 
parts (a) of figures 5 to 10 and by observing the dynamic 
response factors u in table III. The dynamic overstress 
attributable to structural flexibility is the dynamic compo- 
nent of stress in parts (a) of figures 5 to 10 and is represented 
in the response factor u by the amount that this factor differs 
from unit}'. Since response factors as high as 1.97 are ob- 
tained, the dynamic overstress sometimes contributes an 
increment of stress almost as large as the static water-load 
stress. This observation is in agreement with the results of 
other investigators and shows the necessity of using dynamic 
analyses in landing-load investigations. 

The change in wing bending moment due to change in air 
load on the wing is a function of the changes in velocity and 
angle of attack. (See the section entitled “Air-Load Varia- 
tion.”) In these tests the changes in velocity during the 
impacts were small (table III) and the changes in air-load 
bending moment were therefore almost entirely attributable 
to the changes in angle of attack. Since the angle of attack 
is a function of the trim and flight-path angles, large changes 
in angle of attack will occur when the trim and flight-path 
angles change to a great degree. For the relatively small 
changes in trim angle and flight-path angle which occurred 
in these tests (the maximum values being 9.31° and 7.0°, 
respectively) , air-load changes as huge as 0.2 <7 were computed. 
For this airplane, these changes in air-load bending moment 
were of about the same magnitude as the bending-moment 
changes due to the static component of loading (parts (a) of 
figs. 5 to 10) and inclusion of the effects of air-load changes 
in the calculations was therefore necessary. For other air- 
planes with structural and mass characteristics conducive 
to large inertia-load moments, the bending-moment changes 
accompanying a Q.2g change in air load would be small 
relative to the changes in bending moment caused by inertia 
loads. However, for more severe changes in flight-patli 
angle, which should be considered in a design analysis, the 
effects of the change in air load on the bending moment in 
the wing may still be large enough to warrant consideration 
in a design analysis. Further investigation is necessary to 
determine the importance of this air-load variation in design- 
strength calculations. 

In most landing tests the applied forcing functions are 
not easily obtained from the center-of-gravity accelerations 
because of the superposed accelerations caused by structural 
oscillations. But because the fuselage of this airplane is 
connected to the wing by struts located near the nodal point 
of the fundamental mode which represented the greatest 
portion of the wing bending and by a nonstructural fairing 
which neither transmitted nor interfered with the wing oscil- 
lations, the accelerometer in the fuselage was not appreciably 
affected by the wing oscillations in these tests. 


The calculation of the dynamic response of each of the nor- 
mal modes of the seaplane involves solving for the response 
of an equivalent simple spring-mass system to the given 
forcing function. In this calculation it is common practice 
to approximate the forcing function in order to simplify the 
computation. However, the responses used in this report 
were computed from the actual forcing functions because 
trial calculations showed that errors as large as 20 percent 
in the total response of the simple spring-mass system 
could be introduced by use of apparently good simple 
approximations. 

CONCLUSIONS 

Experimental wing bending moments obtained from a 
landing-impact investigation with a small full-scale seaplane 
were compared with analytical results. The effects of the 
variation in air load during the impacts were included in the 
analytical procedure. The responses of the fundamental 
mode were calculated from the recorded time histories of 
the applied hydrodynamic forces. For the seaplane tested 
and the conditions of the impacts encountered, the following 
conclusions may be drawn: 

1. Good agreement between measured and computed 
wing bending moments was obtained when the three com- 
ponents of the wing moment (static water-load moment, 
dynamic water-load moment, and air-load moment) were 
included in the calculations. 

2. The effects of structural flexibility on the wing bending 
moments, represented by the dynamic overstress and under- 
stress, were large, the moment due to the dynamic component 
of the total response being as much as 97 percent of that 
caused by the static water-load component. 

3. Although the changes in seaplane attitude during the 
landing impacts were small, the variation in the air-load 
component of the total moment was of about the same 
magnitude as the static water-load component. Although 
this comparison of changes is not representative of the rela- 
tive importance of the air-load variation in seaplanes with 
structures conducive to large inertia-load components, it 
indicates the probable significance of the effects of air-load 
variation since large changes in seaplane attitude must also 
be considered. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Ya., December 21, 1949. 
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